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Abstract. The Cold Dark Matter paradigm predicts vast numbers of dark matter sub- 
halos to be orbiting in galactic halos. The sub-halos are detectable through the gaps they 
create gaps in stellar streams. The gap-rate is an integral over the density of sub-halos, 
their mass function, velocity distribution and the dynamical age of the stream. The rate 
of visible gap creation is a function of the width of the stream. The available data for 
four streams: the NW stream of M31, the Pal 5 stream, the Orphan Stream and the 
Eastern Banded Structure, are compared to the LCDM predicted relation. We find a 
remarkably good agreement, although there remains much to be done to improve the 
quality of the result. The narrower streams require that there is a total population of 
order 10 5 sub-halos above 10 5 M Q to create the gaps. 



1. Introduction 



LCDM simulations predict about 10 4 sub-halos more massive than 10 6 M Q in a steeply 



rising mass function that constitutes nearly 10% of a typical galaxy's dark halo dSpringel et al 



; iDiemand et al.l 120071 : IStadel et al.l 120091) . The large population of sub-halos is 
not detected in stars, gas, or in any currently detectable annihilation radiation. Logi- 
cally there are three broad possibilities for the deficiency: dark matter sub-halos never 
formed with the LCDM predicted abundance, or, the dark matter sub-halos did form 
and through subsequent evolution were severely depleted, or, the dark matter sub-halos 
do exist but only a few are populated with visible baryons. Any one of these solutions 
has significant implications for galaxy formation theory and cosmology. 

Here we develop a statistical method to count dark matter sub-halos. The gravita- 
tional field of a sub-halo induces orbital deflections in s tellar streams tha t are potentially 

observable as either velocity or positional dis turbances dlbata et al.l2002l : ISiegal-Gaskins & Valluri 



20081 : lYoon et al.11201 ll : ICarlberg et al.ll201 lb . In particular, close passages create gaps 



which are detectable in measurements of the density along a star stream. We provide 
an expression for the average rate at which gaps are created in star streams that are pop- 
ulated with the full LCDM population of dark matter sub-halos. Combining this with 
the predicted length of the gaps gives a predicted relationship between two observable 
quantities, the gap creation rate and the width of a star stream. 



1 



2 



Ray Carlberg 






0.85 0.86 0.87 0. 

fraction of circle 



Figure 1. First panel: The density at the location of closest passage, showing 
decreasing dip depth with increasing impact parameter. Second panel: The depth of 
the density dip as a function of the impact parameter ratioed to the sub-halo scale 
radius for masses increasing from 10 6 to 10 9 M Q . Third panel: The vertical extent of 
the loop for rings at 30, 60 and 90 kpc with the impact approximation prediction for 
60 kpc is shown as the line. Fourth panel: The ratio of the gap length to the vertical 
extent, again for masses increasing from 10 6 to 10 9 M . 



2. The Gap-Rate vs Width Relation 

The rate per unit length at which sub-halos cross the stream at velocity v ± to create gaps 
is, 



r r f „ 

JM Jv ± JO 



(r)N(M)v ± f(v ± )n db dv ± dM, 



(1) 



where b nwx is the biggest impact parameter that can create a visible gap. The various 
factors in Equation Q] are evaluated with numerical si mulations. We t ake th e sub-halo 
mass function N(M) and radial distribution n(r) from ISpringel et al.1 (120081) . We per- 
form a large set of simulations of rings with single halo fly-bys to determine the depth 
and length of the gap that a sub-halo induces. It is quickly apparent that only encounters 
where the sub-halo partially passes through the star stream lead to a noticeable density 
dip. We describe the relation as a ste p function at b ^ qR A M) where R S (M) is the scale 
radius of the sub-halo. We use the ISpringel et al.l d2008h measurement for n(r) (their 
Fig. 11) which finds that the sub-halo density at 100 kpc is 6 times the mean within the 
V50 volume. The simulations provide the necessary relations for the maximum effective 
impact parameter and length of the gap, as shown in Figure [TJ 

The width of the stream is set by the velocity dispersion around a mean guiding 
center motion of the stream. The epicyclic motion of stars around the guiding center 
means that to a good approximation that any gap narrower than the stream width will 
be blurred out by the random motion in the stream. Therefore the stream width is 
the minimum visible gap size and we can eliminate M from our equations through 
w = l(M). Putting all these factors together gives gap creation rate as a function of the 
width of stream, 



Ku(w,r) - 0.059 



n(r)/n Q 



100 kpc 



0.55 



w 



-0.85 



kpc J Gyr 



(2) 



for w in kpc. Equation [2]provides a straightforward testable prediction of the CDM halo 
substructure model, effectively counting sub-halos. 
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Table 1. Gap Statistics and Widths of Star Streams 



Stream 


Gaps 


Length 


Width 


Age 


Rgc 


n/n 


Reference 




# 


[kpc] 


[kpc] 


[Gyr] 


[kpc] 






M31-NW 


12 


200 


5.0 


10 


100 


6 


Carlbera etal. (2011) 


Pal 5 


5 


6.5 


0.11 


7 


19 


24 


Grillmair & Dionatos (2006) 


EBS 


8 


4.7 


0.17 


7 


15 


30 


Grillmair (201 n 


Orphan 


2 


30 


1.0 


3.9 


30 


20 


Newbera et al. (2010) 




0.1 1 10 

width [kpc] 



Figure 2. The estimated gap rate vs stream width relation for M31 NW, Pal 5, 
the EBS and the CDM halo prediction. All data have been normalized to 100 kpc. 
Predictions for an arbitrary alternative mass functions, N(M) oc M~ 16 normalized to 
have the same number, 33, halos above 10 9 M Q as the standard distribution, is shown 
with a dotted line. The predicted relation has a width due to the width of the fitted 
relationships, as well as a halo-to-halo variance. 
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3. Comparison of Predictions and Observations 

The observed stream gap rate and width estimates are summarized in Table 1. The 
results of this paper are largely encapsulated in Figure [2l which plots the measured 
and predicted gap rates against the widths of the streams. The data have been scaled 
for plotting to a galacto-centric distance of 100 kpc using n(r)/no and r in Equation [2] 
Comparison of the CDM based prediction of the gap rate-width relation with the pub- 
lished density data for four streams shows generally good agreement within the fairly 
generous measurement errors. The result is a statistical argument that the vast pre- 
dicted population of sub-halos is indeed present in the halos of galaxies like M31 and 
the Milky Way. 

Question: (Michael Rich) Doesn't the existence of wide binaries limit the number of 
sub-halos that can be present in a galactic halo? Reply: Probably not. The distribution 
of wide binaries is said to be an inverse power (Opik's) law (Longhitano & Binggeli 
2010) to the galactic tidal limit at lpc, with no features. The origin of this distribution 
is unclear. Penarrubia etal (2010) have done some simulations of the survivability of 
wide binaries in the very low velocity dispersion environment of a dwarf galaxy, finding 
that wide binaries should be eroded to a roughly inverse square law beyond 10 4 to 10 5 
au (i.e. nearly a parsec) depending on sub-halo fraction. Therefore it appears that the 
effect will cut in about where the distribution is depleted anyway, so at present there is 
no clear problem from wide-binaries. Carr & Sakellariadou considered this issue in a 
more general context and found the very hard potentials of globular clusters to be more 
of a problem. 

Question: (Joss Bland-Hawthorn) Could the stream clumps be due to "bead" instabil- 
ities? Reply: This is a complex and controversial topic. Kupper, Heggie and others 
have pointed out that the first one or two clumps close to a globular cluster progenitor 
are likely the result of epicyclic pileup. Further down the stream our current full n-body 
simulations find no instabilities, but much depends on the dark matter content. 
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